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Abstract. Production, elliptic flow and the nuclear modification factor of charm and bottom 
quarks are studied in central and non-central heavy-ion collisions at RHIC and LHC using 
the partonic transport model Boltzmann Approach of MultiParton Scatterings (BAMPS). 
^ l' Employing an initial heavy quark yield obtained with PYTHIA the full space-time evolution 

D , of charm and bottom quarks in the quark-gluon plasma (QGP) is carried out with BAMPS, 

^ ■ taking also secondary production in the QGP into account. Only elastic collisions of heavy 

quarks with particles from the medium cannot describe the experimentally observed elliptic 
, flow and nuclear modiflcation factor. However, using an improved Debye screening and the 

^ ' running coupling yields a result which is much closer to data. 
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OO ■ 1. Introduction 

. Due to their large mass heavy quarks are produced at an early stage of the collision and can 

cover - depending on their production point - a long distance through the QGP. Interaction 
along this path and subsequent modifications of the heavy quark can reveal valuable information 
about the properties of the medium. 
^ ' The experimentally observed elliptic flow V2 and nuclear modification factor Raa of heavy 

quarks [H [21 [3] hint at an energy loss which is comparable to that of light quarks. This is in 
contradiction with the expectations which one would draw from the "dead cone effect" [H [5]. 
Whether this large energy loss is due to collisional or radiative interactions - or both (or even 
other effects) - is under investigation (see [3] for a recent overview and comparison with data). 

After a short introduction to BAMPS we will briefly address the production of heavy quarks. 
In ISec. 41 our results on the elliptic flow and nuclear modiflcation factor are discussed and 
compared to the experimental data. 

2. Parton cascade BAMPS 

The partonic transport model Boltzmann Approach of MultiParton Scatterings (BAMPS) OE] 
simulates the full space-time evolution of the QGP by solving the Boltzmann equation, 
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for on-shell partons and pQCD interactions. So far gluons (g) and heavy quarks (Q) are 
implemented with the following interactions: gg — )• gg, gg ggg, ggg gg, gg QQ, 
QQ ~^ 99-, 9Q ~^ 9Q and gQ gQ. Details of the model and the employed cross sections can 
be found in tSlEllH]. 



3. Heavy quark production at RHIC and LHC 

We study the production of heavy quarks in two stages: initially during hard parton interactions 
in nucleon-nucleon scatterings and secondary production during the evolution of the QGP. Due 
to their large mass heavy quark production should be describable in the framework of pQCD. 
Indeed, next-to-leading order calculations [9] are in good agreement with experimental data 
from PHENIX [10]. For the present study, however, we use the event generator PYTHIA [11] 
to determine the initial heavy quark distributions, which agree with the PHENIX data as well. 
Nevertheless, these distributions have large uncertainties due to their sensitivity on the parton 
distribution functions in nucleons, the heavy quark masses as well as the renormalization and 
factorization scales (see [8j for a detailed analysis). For the initial distribution of the gluonic 
medium we use three different approaches: the mini-jet model, a color glass condensate inspired 
model and also PYTHIA in combination with the Glauber model. 

In the following we give a brief overview of our results on heavy quark production in the 
QGP. More details concerning this section can be found in [8j. 

Secondary heavy quark production in the QGP is studied within a full BAMPS simulation 
of Au+Au collisions at RHIC. According to our calculations the charm quark production in 
the medium lies between 0.3 and 3.4 charm pairs, depending on the model of the initial gluon 
distribution, the charm mass and whether a, K = 2 factor for higher order corrections of the cross 
section is employed. However, compared to the initial yield these values are of the order of a few 
percent for the most probable scenarios. Therefore, one can conclude that charm production at 
RHIC in the QGP is nearly negligible. 

At LHC, however, the picture looks a bit different: Here the charm production in the QGP 
is a sizeable fraction of the initial yield and is even of the same order for some scenarios (with 
mini-jet initial conditions for gluons with a high energy density). In numbers, between 11 and 
55 charm pairs are produced in the QGP. 

Bottom production in the QGP, however is very small both at RHIC and LHC and can be 
safely neglected. As a consequence, all bottom quarks at these colliders are produced in initial 
hard parton scatterings. 




4. Elliptic how and nuclear modification factor of heavy quarks at RHIC 

Commonly used observables for investigating the coupling of heavy quarks with the medium are 
the elliptic flow 
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{px and Py are the momenta in x and y direction in respect to the reaction plane) and the nuclear 
modification factor 

""^^ Nundm^p/dpTdy ^""^ 

of heavy quarks at mid-rapidity. A large elliptic flow comparable to that of light partons indicates 
a strong coupling to the medium. On the other hand a small Raa is a sign for a large energy 
loss of heavy quarks. Experimental results reveal that both quantities are on the same order as 
the respective values for light particles [H El [3] . 



As we have recently shown [12] elastic scatterings of heavy quarks with the gluonic medium 
using a constant coupling Og = 0.3 and the Debye mass for screening the t channel cannot 
reproduce the experimentally measured elliptic flow. In order to explain the data one would 
need a 40 — 50 times larger cross section than the leading order one. Of course, this K factor is 
too large to represent the contribution of higher order corrections. However, as we will show in 
this section, the discrepancy with the data can be lowered - even on the leading order level - by 
a factor of 10 by taking the running of the coupling into account and by improving the Debye 
screening. The remaining factor of 4 difference could then indeed stem from neglecting higher 
order effects, which, however, must be checked in a future project. 

The following calculations are done analogously to [131 [H] . An effective running coupling is 
obtained from measurements of e^e~ annihilation and non-strange hadronic decays of r leptons 
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with /3o = 11 - I Nf, Nf = Om our case, and L± = ln(±Q7^^) with A = 200 MeV. If asiQ"^) 
is larger than a™*^ = 1.0 it is set to a™^^. 

Since the t channel of the gQ — >■ gQ cross section is divergent it is screened with a mass 
proportional to the Debye mass m£): 

t t — K rrij-i 

The Debye mass is calculated by the common definition = 47r (1 + Nj /Q) as{t) with the 



running coupling. The prefactor k in Eq. (5) is mostly set to 1 in the literature without a clear 
reason. However, one can fix this factor by comparing the dE/dx of the born cross section with 
K to the energy loss within the hard thermal loop approach to k ~ 0.2 \13\ I14j. 

These improvements lead to an enhanced cross section which also increases the elliptic fiow. 



Fig. 1 shows V2 as a function of the transverse momentum pT for the leading order cross section 
without any improvements, with the running coupling, with the corrected Debye screening and 
with both modifications. The elliptic fiow of the latter reproduces the order of magnitude of 
the data, if the cross section is multiplied with K = 4, which is much smaller than the previous 
employed K = 40 — 50 and lies in a region which could account for higher order corrections. 
However, one has to check if these corrections have indeed a similar effect as a constant K factor 
of 4. Therefore, the calculation of the next-to-leading order cross section is planned for the near 
future and will complement 2 o 3 interactions for gluons, which are already implemented in 
BAMPS [6]. The shapes of the theoretical curve and of the data points are, however, slightly 
different. This is probably an effect of hadronization and decay to electrons, which is neglected 



in this calculation. That is, in Fig. 1 V2 and Raa of heavy quarks on the quark level are 
compared to the heavy fiavor electron data. We plan to take these effects into account in future 
investigations. 



On the right hand side in Fig. 1 the Raa of heavy quarks is depicted, which shows for K = 4 
the same magnitude of suppression as the data. Studies on the Raa of gluons in BAMPS are 
presented in [161 IIZ] • 



5. Conclusions 

We have studied the production and space-time evolution of charm and bottom quarks with 
BAMPS. Charm production in the QGP at RHIC is to a good approximation negligible, but at 
LHC it can reach values comparable to the initial yield. Bottom production during the evolution 
of the medium can be neglected at RHIC and LHC. 
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Figure 1. Elliptic flow V2 (left) and nuclear modification factor Raa (right) of heavy quarks 
with pseudo-rapidity \r]\ < 0.35 at the end of the QGP phase for Au+Au collisions at RHIC with 
an impact parameter of 6 = 8.2 fm. For one curve the cross section of gQ gQ is multiplied 
with a K factor. For comparison data of heavy flavor electrons [3] is shown (see also text). 



The energy loss of heavy quarks in the QGP due to elastic collisions with the leading order 
cross section is very small. Nevertheless, if one improves the calculation of the cross section 
by taking the running coupling and a more precise Debye screening into account, the order 
of magnitude of the experimentally measured elliptic flow and nuclear modification factor can 
be reproduced, if one multiplies the cross section with a K factor of 4. In order to solve this 
discrepancy next-to-leading order contributions must be taken into account, which we plan to 
do in the future. 
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